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Abstract. Dusty galaxies with luminosities in excess of 10^^ Lq have been detected 
out to redshifts 2 ~ 1 by the Infrared Space Observatory (ISO), and to higher redshifts 
using millimetre(mm)- and submm-wave cameras on ground-based telescopes. The in- 
tegrated properties of these more distant galaxies are also constrained by measurements 
of the intensity of the submm-wave background radiation. While it is generally unclear 
whether their energy is released by gravitational accretion or by star formation, cir- 
cumstantial evidence favours star formation. Unless these high-redshift galaxies are 
extremely massive, which is not expected from standard models of galaxy evolution, 
this luminosity cannot be sustained for more than a fraction of a Hubble time, and so 
they are undergoing some sort of 'bursting' behaviour. The interpretation and analysis 
of this population is discussed, and the key observations for deriving a robust history 
of their evolution, which is likely to be the history of starburst activity, are highlighted. 



1 The Evolution of Luminous Dusty Galaxies 

From observations of low-redshift dusty galaxies using the IRAS satellite p8[ 
close to the peak of their restframe spectral energy distributions (SEDs), it is 
known that a similar amount of energy in the local Universe is produced by 
stars in dust-enshrouded and dust-free environments. The comoving luminosity 
density of dusty galaxies is also known to evolve strongly, from the slope of the 
faint counts of /i?j4 5 galaxies at 60|J,m which provide information to z ~ 0.2. 
ISO observations at both shorter jl3j and longer wavelengths confirm that 
strong evolution continues to z ~ 1. 

At longer wavelengths, the redshifted emission from very luminous, high- 
redshift dusty galaxies can be detected in the mm and submm wavebands. Inde- 
pendent surveys made using the 450/850-)J.m SCUBA camera at the JCMT 
have determined the counts of high-redshift dusty galaxies. 1.2-mm surveys us- 
ing the MAMBO detector array at the IRAM 30-m telescope |^ have detected a 
similar population of galaxies. In three cases, the detection of CO emission from 
gas located at the position and redshift of a suspected optical identification (at 
z = 1.06, 2.55 and 2.80) ||l^ provides absolute confirmation of the identification. 
Extremely deep VLA radio images of the survey fields can be used to impose 
constraints on the redshifts and SEDs of the detected galaxies 0,|6| . It is likely 
that the detected galaxies are at z ~ 2 — 3, and there are very few plausible 
low-redshift (z < 1) counterparts. The counts and redshift distributions of these 
distant dusty galaxies can be used to constrain models of galaxy evolution at 
high redshifts. 
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In addition to the detection of individual submm-selected galaxies, the back- 
ground radiation intensity in the mm and submm wavebands traces the 
integrated emission from the entire population of dusty high-redshift galaxies. 
The mm/submm-wave background spectrum has the form I^, oc j/^-^'*, and there 
is no clear spectral break down to 2000 \im. The lack of a break to a steeper slope 
supports the idea that high-redshift galaxies with redshifted SEDs peaking at 
about 1000 jim are still contributing to the background intensity, indicating a 
high maximum redshift of the population of about 10. If the shape of the SED is 
assumed not to evolve significantly, then the luminosity density must evolve as 
PL oc for z 1 in order to generate this background spectrum, 

a result which is independent of cosmology. Such a gently declining high-redshift 
luminosity density is naturally consistent with the incomplete redshift distribu- 
tions of galaxies detected in mm/submm-wave surveys. 

Note that any model of the evolution of dusty galaxies must predict the red- 
shift distribution of submm-detected galaxies correctly. Submm- wave surveys are 
very sensitive to high-redshift galaxies j|] , and it is easy to propose models that 
fit both the observed submm- wave counts and background radiation spectrum at 
the expense of a redshift distribution that is biased far too high. The verification 
of predicted redshift distributions is therefore a crucial test of such models. 

In this paper, the forms of evolution of galaxies that were previously derived 
in the context of far-infrared(IR) and submm-wave data are updated to take 
account of the much greater amount of information that has become available, 
particularly from ISO surveys. The results are similar to, but less uncertain 
than, those derived earlier The results have also been updated to include 

the currently favoured cosmological parameters f?o — 0-3i = 0.7 and Hq = 
65kms^^ Mpc~^ are assumed. 

2 Constraining the Evolution of Dusty Galaxies 

A Baseline at Low Redshifts The luminosity function of IRAS galaxies is 
best constrained at 6G)im [|4|; and information about the same population of 
galaxies is also available at lOOjim [Q. These wavelengths are close to the peak 
of the SED of a nearby (z ~ 0) dusty galaxy for any reasonable dust temperature. 
The ratio of the bright counts at 60 and 100 |lm imply a luminosity-averaged 
dust temperature T ~ 35-45 K. 850-)J,m observations of galaxies detected by 
IRAS 0, with a long wavelength baseline to provide an excellent probe of 
the SED, indicate that T = 36 ± 5K and the Rayleigh-Jeans spectral index is 
3.3 ± 0.2. The population of low-redshift dusty galaxies can be divided into rel- 
atively short-lived warm interacting/starbursting galaxies and long-lived cooler 
quiescent galaxies however, the details of this distinction are relatively 

unimportant for studies of high-redshift galaxy evolution. Any low-luminosity, 
low-temperature dusty galaxies missing from existing surveys do not contribute 
significantly to the luminosity density, even at low and moderate redshifts. 

The form of evolution of the baseline low-redshift far-IR luminosity func- 
tion <Po{L) must be dominated by pure luminosity evolution, that is <I>[L,z) ~ 
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<?o [L/ g{z), 0] , to ensure that the submm galaxy counts and the background radia- 
tion intensity are both predicted correctly. Number density evolution is certainly 
also likely to be involved, but must be dominated by luminosity evolution 
The evolution function g{z) is determined by demanding that the background 
radiation intensity, counts and redshift distributions of dusty galaxies are all in 
agreement with observations. These observations are, in order of increasing red- 
shift, the faintest counts and redshift distributions of 60-}a,m galaxies, deep 
90- and 170-)J,m counts from ISO, the spectrum of background radiation from 
COBE, and the faint counts and limited redshift information of distant galaxies 
detected using SCUBA at 450 and 850|J,m and MAMBO at 1.25 mm. 

Several approaches can be taken to investigate the evolution. The simplest 
is to assume a form for <!>, which implicitly includes details of all the physical 
processes taking place in galaxies, but is fitted to the data without investigating 
the processes in detail . This has the advantage of requiring few parameters to 
model the galaxy SED and the form of evolution, in fact fewer than the number of 
constraining pieces of data. Well-motivated additions of greater complexity can 
thus be introduced to the models as the observations improve, without invoking 
parameters too numerous to constrain reliably and uniquely. A more physically 
motivated approach connects the evolving mass function of galaxies to the asso- 
ciated luminosity function using a prescription for both star formation and the 
fueling of active galactic nuclei (AGN) however, care must be taken 

to avoid getting lost in the space of free parameters. Without an additional pop- 
ulation of short-lived, very luminous galaxies, standard semi-analytical models, 
which include star formation in the gas that cools in galaxy disks, fail to account 
for the observed surface density of SCUBA and MAMBO sources 0. 



Describing Evolution with a Simple Luminosity Function The first re- 
sults derived using this approach Q followed rapidly behind the first results 
of SCUBA surveys ||2^ . A low-redshift 60-)a.m luminosity function was assumed 
[ p4| , an SED was defined by a dust temperature T, and a form of low-redshift 
evolution g(z) = (1 -I- z)'^ was included. T and 7 were determined by requiring 
that the form of the 60-)J.m IRAS counts and the early results of deep 175-|im 
ISO surveys were reproduced; T = 38 ± 4K and 7 = 3.9 ± 0.2 were required 

As discussed above, this temperature is consistent with subsequent SCUBA 
measurements of dust temperatures for IRAS galaxies fl^ , while 7 matched the 
value inferred from optical surveys |2l|] , rather than the value of 7 = 3 that 
is often assumed to describe the evolution of galaxies in the far-IR waveband. 
Note that 7 ~ 4.5 is derived from 15-)im ISO surveys, taking into account the 
complex restframe SED of a dusty galaxy between 5 and lOjim ]3C[| . 

It is now possible to use the much more extensive data from 90- and 175-)J.m 
ISO counts jl^ , and a more popular non-zero- yl cosmology to revisit the results. 
There are no substantial changes; formally T = 37±3K and 7 — 4.05 ±0.15 are 
the latest results, if a Rayleigh- Jeans spectral index of 3.5 is assumed. 

The form of evolution at higher redshifts, too distant for ISO observations, is 
constrained by the background radiation intensity |14| and the counts of SCUBA 
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galaxies. From 2002, mid-IR SIRTF observations should make a major impact 
in this area. Constraints from the background and SCUBA data are rather de- 
generate, although even the first SCUBA data in 1997 provided the tighter 
constraint: see in which various models of high-redshift evolution were con- 
sidered. The so-called Gaussian model, in which g{z) is represented by a Gaus- 
sian in cosmic epoch, provided the best description of the redshift distribution 
of SCUBA galaxies. The more accurate 175- and 850-|lm counts now available 
provide some additional information, and are useful for updating the results. 
Progress has also been made in developing a more appropriate form of g{z), 
which is fully compatible with models of cosmic chemical evolution, and natu- 
rally includes a peak in the evolution function |p7| , p2[ : 

g{z) = {l + z)3/2sech^ [h ln(l + z) - c] cosh^c. (1) 

At low redshifts 7 ~ (3/2) + 2b\/ 1 — sech^c. Using all available observational 
data, the results h = 2.2 ± 0.1 and c = 1.84 ± 0.1 are obtained; see the solid 
line in Fig. 1. These results are similar to those derived earlier but take 
account of revised cosmological parameters and tighter error bars on some of 
the constraining data. The model described above is consistent with all the 
observed background radiation, counts and redshift distributions of galaxies at 
wavelengths longer than 60 )J.m. The dominant source of energy in the Universe 
remains the restframe far-IR radiation of starlight and AGN emission reprocessed 
by dust. 



Describing Evolution Using a Model of Merging Galaxies In an alter- 
native investigation, we took a simple form of the evolution of the merger rate 
of dark- matter halos J0,0| , which adequately reproduces the results of recent 
TV-body simulations |18|. We assumed that a certain redshift-dependent frac- 
tion x{z) of the total mass of dark and baryonic matter involved in mergers is 
converted into energy by nucleosynthesis in high-mass stars with an efficiency 
0.007c2. Note that the same formalism is appropriate for describing the evolu- 
tion of AGN fueling events at the epochs of mergers . The form of evolution 
and normalization of x{z) can be determined by a joint comparison with the 
background radiation intensity and low-redshift IRAS counts. Using an appro- 
priate form of x{z) = giz)/{l + zfl'^ 011, ^ = 1-95 ± 0.1, c = 1.6 ± 0.1 and 
a:;(0) = 1.35 x 10~^ are required in the standard cosmology, assuming the galaxy 
SED discussed above. The resulting history of galaxy evolution is shown by the 
thick dashed line in Fig. 1. 

In addition, the counts of luminous galaxies associated with merger-induced 
bursts of activity can be determined if a fraction F of mergers are assumed 
to generate luminous bursts of duration a. The product _Fcr, which could de- 
pend on redshift, is the function that controls the results. Using information 
derived from the counts of both low- and high-redshift galaxies, the form of 
Fa{z) required to account for the observations can be determined. If the form 
Fa{z) = Fcr(O) exp(az -t- bz"^ + cz^) is chosen, then values of Fa{0) = 2.4 Gyr, 
a = —4.14, b = —0.56 and c = 0.46 are required. 
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Faint Radio and 15-)J.m Mid-IR Counts Both the simple and hierarchi- 
cal models account for all the current observations in these wavebands. If the 
standard form of the far-IR-radio correlation is assumed, with a radio spectral 
index of —0.65, then the 8.4-GHz counts brighter than 10|a,Jy predicted in the 
two models are 1.05 and 0.98 arcmin"^, with slopes of —1.4 and —1.3 respec- 
tively, matching the observed count iV(> S) = (1.01 ± 0.14)(S'/10|lJy)-i-2^^°-2 
[ p3| . The faintest 1.4-GHz counts ^ are also reproduced accurately. If the mid- 
IR SED is described by a power-law f^, oc ly" with a = —1.95 at wavelengths 
shorter than the peak of the SED, then the normalization and general features 
of the deep counts of galaxies determined using 750 at 15|im flS) are also repro- 
duced, including the marked change of slope at flux densities between 0.5 and 
1 mJy. The presence or absence of a PAH emission feature in the SED has little 
effect on the results. In the hierarchical model, the slope of the predicted 15-|a,m 
counts at flux densities between 1 and lOmJy is steeper as compared with that 
in the simple model, in better agreement with the observations. 



3 Conclusions 

These models, which involve a minimum number of free parameters, provide 
a reasonable description of all the counts and redshift distributions of dusty 
galaxies at both high and low redshifts. As additional data, especially more 
complete redshift distributions and very deep SIRTF mid-IR counts, become 
available, more details can be incorporated into the models to reveal further 
information about the properties of evolving distant dusty starbursts/AGN. 
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Fig. 1. The history of star formation inferred using both methods discussed in the 
text. The thick solid and dashed lines represent the simple luminosity evolution model 
and the hierarchical model respectively. The thinner solid lines show the approximate 
envelope of 68% uncertainty in the results of the simple model. The data points are 
taken from a variety of sources: references can be found elsewhere @,|o[. The thin and 
thick dotted lines represent the best-fitting results obtained in previous derivations, 
the modified Gaussian and 35-K hierarchical models respectively. The absolute 
normalisation of the curves depends on the assumed stellar initial mass function and 
the fraction of the dust-enshrouded luminosity of galaxies that is generated by AGN. 
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